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Abstract. We examine the localization properties of the three-dimensional (3D) Anderson 
Hamiltonian with off-diagonal disorder using the transfer-matrix method (TMM) and finite- 
size scaling (FSS). The nearest- neighbor hopping elements are chosen randomly according to 
tij £ [c— 1/2, c+ 1/2]. We find that the off-diagonal disorder is not strong enough to locahze 
all states in the spectrum in contradistinction to the usual case of diagonal disorder. Thus for 
any off-diagonal disorder, there exist extended states and, consequently, the TMM converges 
very slowly. From the TMM results we compute critical exponents of the metal-insulator 
transitions (MIT), the mobility edge Ec, and study the energy-disorder phase diagram. 
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1 Introduction 

From the scaling hypothesis of localization it is well known that in one (ID) and two 
dimensions (2D) almost all states of a non-interacting disordered quantum system are 
localized Delocalized states occur only for special choices of physical parameters. 
An example of such a model is provided by the Anderson model of localization with 
random hopping |^], 

N N 

H^Y.Us\^){J\ + J2e,\^){^\ • (1) 

The hopping elements tij between sites i and j are restricted to nearest neighbors and 
chosen randomly from the interval [c — w/2, c + 'w/2]. Thus c represents the center and 
w the width of the off-diagonal disorder distribution. The onsite potential energies Ci 
are taken to be randomly distributed in the interval [—W/2, W/2]. We set the energy 
scale by keeping w = 1 fixed and use periodic boundary conditions. 

In this work we will concentrate on the properties of the random hopping model 
in 3D. We study the density of states (DOS). Using TMM ||] together with FSS, we 
determine the mobility edge Ec of the MIT separating extended from localized states. 
We find that the off-diagonal disorder is not strong enough to localize all states in the 
spectrum in contradistinction to the usual case of diagonal disorder. 
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2 Density of States 

We have numericaUy computed the spectrum of the Hamiltonian (|^) for 1000 samples, 
system size N = = 10^ and various values of c. In Fig. |l| we show the resulting 
DOS for strong off-diagonal disorder c = and moderate c = 0.6. For c = the peak 
at energy i5 = in the DOS of the finite system is well pronounced. It corresponds 
to a logarithmic singularity in the infinite system |5| as shown in Fig. ^ (left, inset). 
For larger values of c, the off-diagonal disorder becomes weaker and the DOS shows 
various sub-bands which reflect the spectral structure of the completely ordered model 
for w,W = 0. We remark that in the corresponding 2D model the smallest localization 
length and thus the strongest disorder was found at c « 0.25 
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Fig. 1 Left: DOS for off-diagonal disorder with system size 10^ and c = 0; Inset: logarithmic 
behaviour of the DOS near the bandcenter. Right: DOS for c = 0.6; sharp lines correspond 
to the DOS of the ordered system at w,W = 0, scaled by 1/40 for clarity. The binwidth in 
all cases is 0.01. 



3 Determination of the mobility edge 

The application of the TMM for the present case of random hopping relies on a refor- 
mulation of the Schrodinger equation |^ as 

ijL+l^jfcn+l = ~ ejkn)'ipjkn " ^LV'jfcn-l - tf+lkn'4'] + lkn 

— tfkn^J~lkn - tfk+lni^jk+ln — tfkn'^jk-ln (2) 

with '4' jkn corresponding to the wave function at site (j, /c, n), tj-^^ and t^,^ representing 
the hopping elements from site (j — 1, k, n) or (j, k — l,n) respectively to site (j, fc, n) 
and t^ji-^ the hopping element from {j,k,n — 1) to {j,k,n). This equation may be 
written in matrix from as 
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Fig. 2 Left: Reduced localization length A/A/ as function of energy iJ at c = and W — 
as obtained by TMM with 0.1% accuracy. Note the MIT at « 1.28. Inset: Enlarged 
region close to the MIT, only every second data point indicated by a symbol, the error bar 
shown separately. Right: FSS curve for the TMM data from the left panel. 



with 4'n = (■0iin, ■0iMn, ^2in, V'MA/n) denoting the wave function in the nth 
slice, H± the Hamiltonian in the nth shoe, and the matrix tJl connecting the (n — 
l)th with the nth slice. As usual, the localization length A = l/7min is computed 
from the smallest Lyapunov exponent 7min of the 2M^ eigenvalues exp(±7i) of F = 
limK^ooii'l^TK)^^'^^ . Here, tk = TkTk-i ■ ■ ■T2T1. Assuming the vahdity of one-pa- 
rameter scaling close to the MIT, we expect the reduced localization lengths X{M)/M 
to scale onto a scaling curve \{M)/M = f{£,/M) with scaling parameter ^. 

In Fig. H we show the results up to M = 14 at c = and W — 0. The crossing 
between extended states at E < 1.28 and locaHzed states at £' > 1.28 is clearly visible. 
The mobility edge is thus Ec « 1.28. The inset zooms into the MIT region, where one 
can see an M-dependent shift of the crossing point towards smaller E for increasing 
M. Due to this finite size effect, which is supposed to vanish for M large enough, we 
take non-linear and non-universal corrections to FSS Q into account. The FSS curve 
constructed from these data in Fig. || shows two branches similarly to the FSS curves 
in the Anderson model of localization with pure diagonal disorder However, in 
contradistinction to the model with diagonal disorder, we find that the off-diagonal 
disorder at c = is not strong enough to localize all the states. 

Analysing data for larger values of c, we find increasing localization lengths. Con- 
sequently the mobility edge increases with increasing c. In Fig. ^ we show the results 
up to c = 1. For larger c the transition observed for c = in Fig. || becomes less 
easy to detect and we obtain reasonably accurate crossings for large M only. The 
crossing becomes sharp again when varying W instead of E. This allows us to deter- 
mine the critical exponents i>e and vw of the scaling parameter ^ cx |£' — -Ec^^ and 
5 oc |VF — WcY"^ ■ In agreement with recent results in the Anderson model with pure 
diagonal disorder |Q we find ve = 1-61 ± 0.07 and vw — 1.54 ± 0.03. 



510 



Ann. Phys. (Leipzig) 8 (1999) 5 



LU 4 




extended 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Fig. 3 Phase diagram of the Anderson model of locahzation with off-diagonal disorder c 
as obtained from TMM with 1% accuracy. The values of Ec are indicated by diamonds. 
Squares represent the numerically determined upper band edge Eb for a system with 10"^ 
sites averaged over 1000 samples. The straight dashed line is the band edge Emax for the 
infinite system corresponding to a constant t = c+ 1/2. The difference between Eb and -Bmax 
occurs due to exponential tails of the DOS which cannot be resolved by our calculations. 
States with \E\ < Ec are extended (dark gray), whereas states with Ec < \E\ < Eb are 
localized (light gray). 



4 Conclusions 

We have investigated the DOS and the locahzation properties of the Anderson model 
of localization with ofF-diagonal disorder. Our results show that in 3D a finite part 
of the spectrum close to the band center remains extended even for the strongest 
off-diagonal disorder c = 0, whereas in ID and 2D the off-diagonal disorder has non- 
localized states only for E — 0. FSS is possible and we determine the mobility edge 
separating extended and localized states in the off-diagonal disorder vs. energy phase 
diagram. The obtained critical exponents of the MITs are in agreement with prior 
results for the model with only diagonal disorder ^ and convincingly support the 
universality. 
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